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We theoretically investigate a photonic Kagome lattice which can be realized in microwave cavity 
arrays using current technology. The Kagome lattice exhibits an exotic band structure with three 
bands one of which can be made completely flat. The presence of artificial gauge fields allows to 
emulate topological phases and induce chiral edge modes which can coexist inside the energy gap 
with the flat band that is topologically trivial. By tuning the artificial fluxes or in the presence 
of disorder, the flat band can also acquire a bandwidth in energy allowing the coexistence between 
chiral edge modes and bulk extended states; in this case the chiral modes become fragile towards 
scattering into the bulk. The photonic system then exhibits equivalents of both a quantum Hall 
effect without Landau levels, and an anomalous Hall effect characterized by a non-quantized Chern 
number. We discuss experimental observables such as local density of states and edge currents. We 
show how synthetic uniform magnetic fields can be engineered, which allows an experimental probe of 
Landau levels in the photonic Kagome lattice. We then draw on semiclassical Boltzmann equations 
for transport to devise a method to measure Berry's phases around loops in the Brillouin zone. 
The method is based solely on wavepacket interference and can be used to determine band Chern 
numbers or the photonic equivalent of the anomalous Hall response. We demonstrate the robustness 
of these measurements towards on-site and gauge-field disorder. We also show the stability of the 
anomalous quantum Hall phase for nonlinear cavities and for (artificial) atom-photon interactions. 



I. INTRODUCTION 

Phenomena associated with flat bands in fermionic and 
bosonic systems have drawn attention over the last three 
decades [THE]. The best known example is perhaps the 
fractional quantum Hall effect [9]. If a single-particle 
band is dispersionless in one direction, electrons are lo- 
calized in that direction causing an element of the effec- 
tive mass tensor to diverge. If the band is completely 
flat in the entire fc-space, then heavy degeneracy appears 
and the density of states diverges. Such singularities in 
the density of states are expected to produce anomalous 
behaviors in physical properties including transport phe- 
nomena and optical response. If a wavepacket is created 
in the flat band system, the group velocity will automat- 
ically vanish and very strong backscattering prevents the 
packet from moving. In fact, this scenario takes place 
on the Kagome lattice where localized states on hexagon 
rings occur as a result of the alternating sign of the wave- 
function amplitude jTOj [11] . 

The topological properties of the band structure are 
now essential in approaching the problem of Chern insu- 
lators 12J or topological insulators [USE]. Bandstruc- 
ture topology has been first discussed in the context of 
the integer quantum hall effect by Thouless et al. |15j . 

More precisely, the set of energy eigenstates that form 
an isolated band is described by the first Chern number 
|16j . a topological invariant associated with the band. 
The Chern number has direct physical consequences: for 
example, a completely filled (electronic) band has a quan- 
tized Hall conductance, corresponding to the existence 
of chiral edge modes [17]. It is possible to view this in- 
variant as the flux of the Berry curvature [18) through 



the first Brillouin zone [15] . The Berry curvature and 
the Chern number have become important mathematical 
concepts to classify topological phases in connection with 
band structures. The Chern number has previously been 
accessed through the conductance and chiral edge states 
of quantum Hall systems, and the total Berry phase as- 
sociated with a single Dirac point has been inferred from 
transport in graphene [20]. Measuring the local Berry 
curvature constitutes an actual challenge in experiments. 

In this paper, we investigate anomalous Hall phases of 
light on a specific Kagome lattice with artificial magnetic 
fluxes, which has been introduced in the context of pho- 
tonic lattices [UJ ; see Fig. [I] The effect of time-reversal 
symmetry breaking in connection with this Kagome pho- 
ton lattice has been previously studied at the level of a 
few site geometry [21] . We show below that this lattice 
with artificial gauge fields allows for the existence of an 
anomalous quantum Hall phase, or a quantum Hall phase 
without Landau levels |12j . 

Photonic lattices based on arrays of circuit quantum 
electrodynamics (cQED) superconducting elements [52] 
have been recently realized experimentally [531 HI] • Such 
photon cavity lattices are predicted to exhibit interesting 
many-body phenomena |25[ 126] including a superfluid- 
Mott transition of light |2"71 - I35] . Bose-Hubbard models 
with attractive interactions [35] , fractional quantum Hall 
physics [37-39! and interesting dynamics 10 in partic- 
ular in ring geometries [41] . Recently a Rabi model has 
been shown to present a Z 2 Ising universality class quan- 
tum phase transition between two gapped phases, which 
are not Mott insulator phases [42] . An anomalous quan- 
tum Hall phase for light accompanied by unidirectional 
photonic chiral edge states has been predicted in photonic 
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FIG. 1. (Color online) The photonic Kagome lattice with 
three sites per unit cell A, B,C. Artificial fluxes can be 
threaded through the triangular plaquettes |21j. Top left: lat- 
tice periodic in the x direction, with "line" boundaries com- 
posed of A, B sites only. The blue parallelogram encloses a 
one-dimensional superunit cell which generates the entire lat- 
tice. Top right: a projection of the Brillouin zone onto the 
k x = k\ direction. Bottom left: at </> — n/6 the system ex- 
hibits a flat, topologically trivial middle band, while the lower 
and upper bands are characterized by finite Chern numbers. 
Bottom right: the middle band becomes dispersive if one de- 
viates from it/6, e.g. for (f> — n/4. Here, we have considered 
a cylinder geometry which allows for two chiral edge modes: 
right-moving at the lower edge (purple), left-moving at the 
upper edge (orange); corresponding chiral currents depicted 
on the lattice. As discussed in Sees. IV B and C, the chiral 
edge modes can also be detected in the square geometry [23] . 



crystals 43, 44] and confirmed experimentally [15]. In ad- 
dition, two-dimensional photonic equivalents of topolog- 
ical insulators have been envisaged (46J [47] . Topological 
properties of optical systems might help to implement 
robust photonic devices [15] and to realize invisibility 
cloaks [48 . Artificial gauge fields in photonic systems 
have attracted growing attention [49, 50J . 

Our primary goal is to investigate the effect of the 
topologically trivial flat band (with a zero Chern number) 
on the robustness of photonic chiral edge modes. In our 
realization of the Kagome lattice, neighboring hexagons 
are subject to the same magnetic flux. The magnetic flux 
opens a gap at the Dirac points and breaks time-reversal 
symmetry allowing to stabilize topological phases char- 
acterized by two bands acquiring non-zero Chern num- 
bers (see Fig. [I]). The middle band can be made flat 
or acquire a finite energy dispersion for certain values 
of the artificial gauge fields (see Fig. [T]) or if disorder 
is present; therefore, edge modes can coexist with ex- 
tended bulk states at the same energy. In this situation 
edge modes are fragile towards leaking into the bulk re- 
sulting in an anomalous Hall phase with a non-quantized 
Chern number [51 . A similar effect has occurred in a 



toy model on the honeycomb lattice [55]. We investigate 
the evolution of unidirectional edge modes of light in the 
presence of a localized scatterer by tuning the artificial 
gauge fields and by introducing disorder. We present a 
manner to determine Berry's phases of photon wavepack- 
ets through an interference experiment. There exists a 
previous proposal to access the Berry curvature through 
a measurement of group velocities of wavepackets in opti- 
cal lattices [53| . We discuss a direct measurement of the 
Chern number of a Bloch band, which does not require 
summation of contributions across the Brillouin zone. 
Underlying these methods is the implementation of uni- 
form synthetic gauge fields through tunable cavities with 
time-dependent frequencies, and the realization of Lan- 
dau levels in photon systems. It shall be noticed that a 
similar scheme has been recently realized in cold atom 
systems [54] and that cavities with tunable resonances 
have been already implemented in cQED (55] [56] . 

In addition, the quest for exotic phases on the Kagome 
lattice represents an active subject of research [U I57r[6"6"] . 
Topological phases and their supporting artificial gauge 
fields are potentially realizable in cold atom systems |67F 
[72]. The Hofstadter spectrum [73 J on the related dice 
lattice has been realized in GaAlAs/GaAs systems |74) . 
The Hofstadter spectrum for the Kagome lattice has also 
been studied with superconducting wire networks |75) . 
Finally, the Kagome structure naturally appears in real 
materials [75J [77] and recently it has been implemented 
in cold atom systems [75]. 

The paper is organized as follows. We introduce the 
tight-binding model of a Kagome lattice in the presence 
of artificial gauge fields in Sec. [n] Explicit solutions for 
the edge mode wavefunctions in cylinder geometry and a 
brief discussion of the bulk-edge correspondence are the 



subjects of Sec. Ill Following this, in Sec. IV we present 



observables specific to the anomalous Hall phase, namely 
the bulk polarization, the edge currents, and the local 
density of states. Sec. [V] further details on a method to 
measure Berry's phases of wavepackets in both clean and 
disordered systems. The method requires a uniform syn- 
thetic magnetic field, on whose realization in a photonic 
lattice we elaborate. Additionally, we show a method 
to directly find the topological Chern number of a given 
Bloch band by counting energy levels of the system in an 
artificial magnetic field. In Sec. |VI[ we discuss in detail 
the anomalous Hall effects of light in relation to Chern 
numbers. In the particular case of disordered systems, 
we present a real-space calculation of the Chern number 
which will reveal a close analogy between tuning the ar- 
tificial gauge fields and the presence of disorder in the 
system (in cQED lattices, disorder can appear either as 
an on-site scalar or as a vector potential). In Sec. VII 
we discuss the role of symmetries and the local stability 
of Dirac points in the three band system. In Sec. |VIII| 
we study the effect of interactions in QED cavities 79 , 
such as Bose-Hubbard [80] and Jaynes-Cummings [51] . 
on the topological phase. The Appendices are dedicated 
to technical details. 
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II. TIME-REVERSAL SYMMETRY BREAKING 

Below, we introduce the Kagome lattice with artificial 
gauge fields. An effective tight-binding Hamiltonian for 
photons with the possibility to break time-reversal sym- 
metry can be realized in a cQED system by coupling su- 
perconducting waveguides to nano-Josephson circulators 
(rings) 21j. The Kagome lattice is formed by triangular 
and hexagonal plaquettes arranged as shown in Fig. [T] 
The resulting structure has three sites per unit cell, which 
we denote by A,B,C. Letting a be twice the length of 
a bond, we pick the following lattice vectors A 12 and 
reciprocal lattice vectors gi 2 : 



A 1 = a (l,0), A 2 



47T 

\/3a 



(1,0), g 2 




(1) 



The Brillouin zone, see Fig. [T] is hexagonal and contains 
the common points of high symmetry T = (0, 0) and 

We will use the following lattice coordinates: the di- 
rections along Ai, A 2 in real-space corresponding to 
gi , g 2 in momentum space (top right panel of Fig. [IJ , 
where hats denote unit vectors. In this notation, Ai 
coincides with the x direction, x. We denote sites on 
the lattice by a pair of integers m = (mi,m 2 ), and co- 
ordinates on the lattice in terms of lattice vectors by 
r m = TOiAi +m 2 A 2 . Continuous coordinates are repre- 
sented by a pair of real numbers r = (ri, r 2 ) in terms of 
unit vectors, r = i\ Ai + r 2 A 2 ; momenta will be written 
as k = fcigi + fc 2 g 2 . 

In a cQED system, a given waveguide of length in 
the millimeter range typically supports a single photonic 
mode. Furthermore, the photons which travel in these 
one-dimensional waveguides do not carry a polarization 
label and must be thought of as simply excitations of mi- 
crowave resonators [82 . The hopping of photons from 
one waveguide to another results in a translationally in- 
variant tight-binding Hamiltonian. Following Ref. |21j . 
a honeycomb array of waveguides with nano-Josephson 
circulators can be equivalently reformulated as the pho- 
tonic Kagome lattice of Fig. [T] The Hamiltonian can be 
Fourier transformed, and we obtain: 



(2) 



keBZ 



and the spinor tp\^ = ^a^ k , o(gk' a Ckj contains the cre- 
ation operators on each sublattice, and 



Huj 2|<|e^cosai 
-74 : : | 2|i|e-^cosai hw 

2|i|e l<?i cosa 2 2|i|e"^ cos a 12 



2\t\e~ 1 ^ cosa 2 
2|t|e^cosai2 



We have defined three dimensionless functions of momen- 
tum, 

Ai - A 2 



(3) 



Ai A 2 

«i(k) = k. — , a 2 (k) = k. — , «i 2 (k) ee k 

(4) 

In the effective tight-binding model particles hop be- 
tween nearest-neighbor sites with a complex hopping in- 
tegral |*|e ±i '*. Photons acquire a phase 30 around a trian- 
gular plaquette and a phase of —6(f) around the hexagonal 
plaquette, amounting to zero total flux in the parallelo- 
gram unit cell of area |Ai x A 2 |, which is characteristic 
of the anomalous quantum Hall effect without Landau 
levels as introduced by Haldane [12!. The particles are 
also subject to an artificial on-site potential induced by 
the (electromagnetic superconducting resonator) waveg- 
uide frequency u, and an artifical gauge field giving rise 
to complex hopping integrals and breaking time-reversal 
symmetry. It shall be noted that both the effective on- 
site potential Htu and the hopping strength \t\ should 
be obtained rigorously after integrating out coupling el- 
ement degrees of freedom (in our case, a Josephson ring 
coupling three resonators) [ST] . In a typical cQED exper- 
iment, the energy of an incoming photon is much larger 
than the size of the energy gap between bands: \t\ is ex- 
pected to lie below 100 MHz and the frequency uj is in 
the GHz, or microwave, range. 

The rotational symmetry of the lattice is not strictly 
necessary and anisotropies are possible; their effect on 
the bandstructure is discussed in Sec. IVHI In the rest of 
the paper we shall focus on the isotropic case. 

The Hamiltonian of Eq. [3] is a band Hamiltonian obey- 
ing Bloch's theorem. To determine the Bloch energies 
E(k) for the three bands let us denote U(k) = —(E(k) — 
Hu))/\t\, and a(k) = cosai(k) cosa 2 (k) cos(«i 2 (k)). We 
obtain the following equation for the eigenvalues 

U 3 (k) - 4Z7(k)(2a(k) + 1) + 16cos(30)a(k) = 0. (5) 

For special values of the phase 0, destructive interference 
confines the wavefunction to the hexagonal plaquettes, 
which renders the respective band completely dispersion- 
less. We obtain such a solution with energy U(k) = 
2 cos 30 whenever cos(30)(cos30 — 1) (cos 30 +1) = 0. 
The middle band is flat whenever = | + where m 
is an integer. When <f> = ^y 1 , the upper or lower band is 
flat, touching the middle band at the T point; the other 
two bands form Dirac cones at the six corners of the Bril- 
louin zone. For example, the bandstructure for the cases 
= ^ and f (which will be thoroughly discussed in the 
text) is depicted in Fig. [T] along with edge modes arising 
from a finite geometry to be discussed below. The edge 
modes appear because by tuning the phase <fi one breaks 
time-reversal symmetry at the Dirac points. The band- 
structure of Eq. ([5]) is invariant to a shift of ^ in 0; we 
then restrict to € [0, ^p] in the further discussion. 

Time-reversal symmetry is kept if the phases accumu- 
lated along any closed loop on the lattice amount to an 
integer multiple of tt |21j . This condition is equivalent 
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FIG. 2. (Color online) Two-dimensional phase diagram of a 
particle (photon) with energy E on the Kagome lattice defined 
with hopping phase <f>: by tuning (E, 0) the photon system 
is either characterized by extended bulk states ( "bulk band" ) 
or by a "band gap" where the system exhibits unidirectional 
edge modes but no bulk states. The F points represent states 
in the flat band (the blue points at E — hw — correspond 
to a middle flat band separated by finite energy gaps from 
the lower and upper bands; black points at E — hu) = ±2|t| 
correspond to flat upper or lower band). The D points are 
degeneracy points between pairs of bands at the Dirac points 
K . The dashed lines are explained in the main text. 



to requiring that J^Lk and its time-reversed counterpart, 
^k*, are the same up to a gauge transformation on the 
spinor -0k in Eq- In our case, if 30 is an integer mul- 
tiple of 7T, then the Hamiltonian is time-reversal symmet- 
ric. The bandstructure properties can be summarized in 
the (E, 4>) phase diagram of Fig. [2j The phase diagram 
represented is periodic in 0. A particle created at energies 
in the "band gap" regions cannot be in a two-dimensional 
Bloch state of the Hamiltonian. States exist in this en- 
ergy interval only in a finite geometry with boundaries, as 
we show in detail in Sec. Mil Such a state would be con- 
fined to the edges of the sample and propagate chirally 
(chiral edge mode). The "bulk band" phase is embodied 
by bulk extended states. When photonic bulk extended 
states and chiral edge modes coexist at the same energy 
E, the edge excitations may scatter into the bulk; this 
situation coincides with an anomalous Hall phase. The F 
points in the diagram of Fig. [2] correspond to flat bands. 
Dashed lines from the F points show the possible phase 
transitions: horizontal blue line - from the flat band to a 
band gap by changing the particle energy; vertical blue 
line: from the flat band to a bulk band by changing the 
phase (j). The D points represent band degeneracies at 
the Dirac points K ± . Moving along the vertical red line 
(a change in 0) lifts the degeneracy and a finite band gap 
is opened. A particle with a slightly larger energy would 
be in a bulk band state. 

At = 0, the system is gapless: there is a quadratic 
touching with a flat band at the T point, corresponding 
to F in Fig. [2j and Dirac cones at the K points, cor- 
responding to D. We will be focusing on the system in 



the vicinity of = |, at which the system exhibits a 
flat band. Changing the energy of the particle at the F 
points (equivalent to moving along a horizontal line in 
the figure) will take us from the flat band into a band 
gap. Changing at the F point (equivalent to moving 
along the vertical line) will take us from the flat band 
into a bulk-band with finite bandwidth. 



III. CHIRAL EDGE MODES 

In this section we make explicit the correspondence 
between a bulk quantity, the Chern number, and the chi- 
ral edge modes at the edges of a sample. It should be 
noted that even though we choose the cylindrical geome- 
try below for mathematical convenience, the edge states 
can be detected in a square geometry; this point will be 
illustrated below in Figs. 5 and 6. 

Given a translationally invariant two-dimensional lat- 
tice Hamiltonian whose eigenstates arc the Bloch states, 
the existence of chiral edge modes is signaled by the non- 
zero values of the Chern number i/W corresponding to 
the Bloch band |nk) [15], 

y{n) = iX/ 2k ( 5kX ^ ( " )(k) ) ' (6) 

where the vector field (k) is the Berry gauge poten- 
tial associated to the n th Bloch band, 

«W(k) = -i{nk|at|nk). (7) 

The integral in Eq. (|6| is over the entire surface of the 
Brillouin zone, so it is a summation over all of the states 
enclosed in the n th band. 

When the bands are separated by an energy gap 
throughout the Brillouin zone, the band Chern num- 
bers of the lower, middle, and upper bands, respectively, 
can be shown numerically to be — sgn (sin 30) , and 
+sgn (sin 30), where 6 [P>ir]' ^ ne sm (30) function 
was chosen to obey the 27r/3 periodicity of the bandstruc- 
ture with respect to the phase 0. In particular, this is 
consistent with a recent result that implies that the flat 
band of the Kagome lattice, if isolated from the other 
bands by energy gaps (as is the case for = ^), will be 
non-topological [53"] . 

The Chern number of Eq. ^ is a topological index 
whose value can only change if by variation of parame- 
ters the bands touch at some degeneracy points in the 
Brillouin zone and subsequently reopen the energy gap 
[15] , As the parameter is varied, for values 0, |, 
two bands touch at all corners of the Brillouin zone form- 
ing Dirac points; the remaining pair has a degeneracy at 
the r point. At the transition upper and lower bands 
exchange Chern numbers leaving the middle band topo- 
logically trivial throughout. The stability of Dirac points 
and the nature of the exchange of Chern numbers is ex- 
plained in Sec. |VII| We now make the correspondence 
between the bulk quantity and the edge modes explicit. 



5 




FIG. 3. Effect of different boundary conditions on the disper- 
sion of edge modes in the proximity of the dispersive middle 
band for a system with (j> — J: left, two "line" boundaries; 
middle, one "line" edge, the other "armchair"; right, both 
boundaries are "armchair". The bottom panels exhibit the 
minimal ladder lattices that exhibit the three boundary con- 
ditions. All three ladders shown support edge modes. 

Consider an analogous fermion system at zero temper- 
ature: all states with energies lower than the chemical 
potential \x are occupied. If [i is set to an energy gap, the 
sum of band Chern numbers for the occupied bands is an 
integer quantity, possibly zero, which is proportional to 
the Hall conductivity [15] . The principle of bulk-edge cor- 
respondence states that this sum counts the chiral edge 
modes [17] supported at the given chemical potential. 
Since the argument is independent of statistics, it applies 
equally well for photons; we will then replace the notion 
of "chemical potential" by "energy of injected particle" . 

Let us place the Kagome lattice in a cylinder geom- 
etry: a particle injected at an energy inside the energy 
gap can belong to one of two states localized at the edges. 
We now find an analytical solution of the edge wavefunc- 
tions. We consider a cylinder with line boundaries along 
the periodic Ai direction (boundaries containing A, B 
sites only - see Fig. [3]). The cylinder is generated by 
translating the superunit cell (blue parallelogram in Fig. 
[I]) along the lattice vector Ai. The superunit cell is ob- 



the first three of which hold for to 2 spanning the supe- 
runit cell "bulk" between and m2 im „, and the other 
three are equations at the boundary. There exist solu- 



tained by translating the triangular plaquette along A2. 
We number plaquettes along the superunit cell by integer 
rn 2 = 0, 1, ...,m 2max , as in Fig. [T] We Fourier transform 
the operators d m =(mi,m 2 ) along the periodic k\ direction 

* mi 

where a denotes the sublattice A, B or C; N\ is the num- 
ber of repeated superunit cells making up the cylinder; 
and integer m,2 indicates the triangular plaquette along 
the superunit cell, as shown in Fig. [T] The problem has 
become to diagonalize the one-dimensional Hamiltonian 
describing a superunit cell: 

,ye kl = \t\ 2e^4 m2 a Am2 cos (fc^) 

m 2 

+e^ fel H m2 a Bm2 + e^ k ^al m2 a B , m2+1 

+e l ^ k ^a\ m2 ac m2 +e^-^ia\ m2 a c , m2 „ 1 + h.c. 

(9) 

This form reduces to the one in Eq. ^ when Fourier 
transformed around the remaining direction fc 2 . To ob- 
tain the edge states, let us take the general form of an 
eigenstate 

a=A,B,C m 2 

where |am 2 ) is a ket localized on sublattice a of the 
m| triangular plaquette along the super-unit cell. The 
Schrodinger equation 

Jft 1 |*(fci)) = S*(*i)l*(*i)> (11) 

yields the following linear system for the components of 
the wavefunction ip am2 = "0am 2 (^i) and the energy dis- 
persion U = U(ki) = - (E(h) - huj) /|i|, 



U^Bm 2 


(12) 


U^Am 2 


(13) 


U 1pC,m 2 


(14) 




(15) 




(16) 


U i>Am 2m!tx , 


(17) 



I 

tions which are exponentially suppressed with distance 
from the boundary of the sample, which we put in the 



e^i+^Vc,m 2 -i + 2 cos (fcr|) e^VUm, + e" ifel t+^ Cm2 = - 
e-^S-'fycm,-! + 2 cos (fci |) e l H Bm2 + e^^Uc m2 = - 
e *f-^ Sm2 +e- ifcl t+^ Am2 + e - ife ^-^B, m2 +i + e*t+^A,m 2 +i = - 

e ik ^ +i ^ c ,m 2 ^ + 2 cos e-*+1> Amtm = - 

e-^-^cm^ + 2 cos (ferj) e 4 ^B™ 2max = - 



G 



form (a denotes sublattice) 



A + 2 , A_ 2 



(18) 



The condition that the two edges of the cylinder are line- 
shaped amounts to requiring 



ipc,-i = Tpc,?n 2n 



0. 



(19) 



A pair of chiral edge modes localized at the boundaries 
of the sample arises: 



' k A cos 
E. = 2\l <ch ( J and A + = - 



V 4 2 J 



E_ 



-2|t| cos 



kia 



sin 



and A_ 



2 



(20) 



If |A±| < (>)1, then the state is confined to the top 
(bottom) boundary. The states identified to be at the top 
(bottom) travel in the +(— )Ai direction, so "chirally". 
The edge modes along with the bulk states for a flat 
band system with <j) — ^ and for a dispersive middle 
band system at <j) — § ar e plotted in the bottom panels 
of Fig. [T] The simple boundary condition chosen gives 
an energy dispersion similar to that of a one-dimensional 
tight-binding chain of lattice constant |. The factors 
A± describing the spatial suppression of the wavefunction 
into the sample bulk can become singular or 0, in which 
case the normalized wavefunction has weight inside of 
the bulk and is pinned to one of the two edges (see Figs. 
[I]and[4]). For other boundary conditions (see Fig. [3] for 
the equivalent "armchair" condition), one can arrange 
that the edge mode and the middle bulk band overlap 
in energy, but not in momentum, or in both energy and 
momentum. States overlapping in energy can scatter into 
each other: in the "line" boundary conditions, a spatially 
localized, <5-function, impurity can scatter the edge mode 
into the bulk, as we show in Sec. IV C 

For the case of a Kagome ladder, the counter- 
propagating edge modes exist as long as there are two 
distinct edges, of either the "armchair" or the "line" type. 
Therefore, any ladder will need to have at least one row 
of hexagonal plaquettes. The configuration with the least 
sites that sustains edge modes is shown in the bottom left 
of Fig. [3] This specific ladder with two "line" boundaries 
has the property that the flux per unit cell is 0. 

An intriguing property related to the bulk-edge corre- 
spondence is the phenomenon of polarization in a topo- 
logically non-trivial system. We develop on this in the 
context of a photonic system in Sec. |IV A| 



IV. OBSERVABLES 

In this section, we elaborate on observable quantities 
in the cQED based photon lattice exhibiting two differ- 




FIG. 4. The quantities A+ (solid line), and XL (dashed 
line) defined in Eqs. ( 20 1 describe how the wavefunction of 



the edge mode decays into the bulk. The reference line is at 
X± = 1. For |A±| > 1 the respective branch E±(ki) is located 
at the bottom edge. Conversely, if |A±| < 1, it is located at 
the top edge. This results in the color scheme for the edge 
modes on Fig. [I] 



the polarization of a topologically nontrivial bulk in the 
context of photonic systems and draw connections to the 
local density of states in Sec. |IVB| Sec. |IV C| is dedicated 
to the lattice current density. Let us emphasize that the 
local density of states and lattice current densities are in 
principle accessible in square geometries. 



A. Polarization 

The phenomenon of polarization is a signature of the 
fact that it is impossible to choose a smooth gauge for 
the Bloch wavefunctions across the Brillouin zone in a 
topologically non-trivial system. Polarization has been 
introduced in the context of electronic Chern insulators 
in Ref. [S3], and it is related to the Hall conductivity. 
In a photonic or cold atom system, polarization may be 
observed in the dynamics of wavepackets, as we are about 
to describe. 

Let us consider the one-dimensional system described 
by the Hamiltonian of Eq. This time we close 
the cylinder around both directions Ai and A2 into 
a torus, such that J4fk ± now describes a periodic one- 
dimensional lattice along the A2 direction, in terms of 
the periodic parameter k\. It is possible to find for this 
one-dimensional periodic system a set of Bloch states 
|nk) = |n fc 2 ;fci) with a smooth choice of gauge every- 
where as a function of fc 2 and for each value of k% . This is 
equivalent to the possibility to construct Wannier states 
which are exponentially localized along the A2 direction, 
which is a classic result of Kohn [55] . For now let us focus 
on a single band |nk) = Ifcife). The maximally localized 
Wannier states take the following form |86| : 

|W(fc ls iZa))= (21) 
* 'No 



ent Hall phases for light In Sec. IV A we discuss where N2 is the number of triangular plaquettes in the 
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quasi one-dimensional system; R 2 is the real coordinate 
of the center of the Wannier function along the A2 di- 
rection, which is made precise below; &i f2 (ki,k 2 ) are 
the components along the directions Ai,A2 of the Berry 
gauge field introduced in Eq. and 



versus localized character of bulk states can be tested by 
measuring the local density of states. Because phase in- 
formation is absent from the density, it is not possible to 
measure the polarization of Sec. |IV A| Given a complete 
set of states |^), each of which is an eigenstate of energy 
Ei& of the tight binding Hamiltonian of Eq. ([2| , 



1(h) 



V3a 



dK 2 M 2 (h,K 2 ). 



(22) 



(26) 



Since the Bloch functions can and have been chosen pe- 
riodic and smooth in the k 2 direction, both of the com- 
ponents of the Berry gauge field ^1,2(^1, k 2 ) are periodic 
functions of k 2 . 

The coordinate of the center of the Wannier functions 



one can define the local density of states in terms of the 
energy E and the two dimensional real coordinate on the 
lattice r 



p(E,r) 



drifts depending on the parameter k\. 
value of the position operator r 2 is 

(W(h,R 2 )\r 2 \W(h,R 2 )) = R 2 



The expectation 



6(E - E 9 )\(r\V)\' 



(27) 



l(h) 
(4tt/V3) ' 



(23) 



This quantity is by definition normalized to the number 
of sites in the system, 



The net displacement of the center of the Wannier state 
when the parameter h winds around its period equals 
the Chern number: 



dE^p(E, 



= 37V, 



(28) 



Ar 2 



1 



(4tt/\/3) 



(4tt/V3) 
' {4n/y/3) k 



(l Uir/VSa 

la / 

dn 2 ( . 



1(0) 



(47r/^ 



dk.M(h,k 2 ) 



3a, n 2 



-va. 



where m = (mx,m 2 ) indexes the sites, r m = miAi + 
TO2 A2, and iV denotes the number of unit cells and there 
are 3 sites per unit cell. We plot the local density of 
states in the same cylinder geometry with "line" edges 
.of Sec. |III| We have approximated the ^-function by a 
'2 (0, K 2 ) JLorentzian 



(24) 



9BZ 



lim — - — 

A E ^I1 7T (E - 



6{E-E*), (29) 



In the last line we have used the fact that 3S\ 2 (hi k 2 ) is 
periodic in k 2 to reexpress the drift of the center of the 
Wannier state as the line integral of the Berry gauge field 
around the boundary of the first Brillouin zone dBZ, 
which is exactly the Chern number v. This condition 
summarizes the fact that a smooth gauge cannot be cho- 
sen for the Bloch functions in a band with non-zero v. If 
this were true, then the Berry gauge field would be pe- 
riodic in both directions, and the Chern number would 
vanish. We arrive at the following boundary condition 
for the Wannier states, 

W (ki + iTT/V3a, r 2 ) ^ = |W(*i,r a - va)). (25) 

In particular, this tells us that if a wavepacket con- 
structed out of such states is accelerated in the Ai direc- 
tion, then its center will drift along the A2 direction. A 
wavepacket constructed within a band of Chern num- 
ber would return to its initial coordinate on the r 2 axis 
upon a complete revolution of k% . We shall come back to 



where the width is set to A.g/|t| = 0.01 and a system of 
23 x 23 sites has been numerically diagonalized. In the flat 
band model with (f> = ~ , as shown explicitly in |21j , there 
exists a highly degenerate band of zero-energy states lo- 
calized in the hexagonal plaquettes due to destructive 
interference effects. When the flux (j> is detuned from |, 
the middle band becomes dispersive and the states leak 
out of the hexagonal plaquettes and become extended 
(for (j> = j, see Fig. 5(c) I. Finally, if disorder is intro- 



this point in Sec. VB 



B. Local density of states 

The localization of chiral modes at the edges of the 
sample, the existence of a flat band, and the extended 



duced into the flat band system at (f> = ^, the states 
inside the middle band become delocalized due to dis- 
order (Fig. |5(d)[ ). Intra-gap edge modes have non-zero 
spectral weight only close to the edges of the sample. 



C. Lattice currents 

A measurement of the local density of states p(E, r) 
cannot detect the chirality of edge modes. The simplest 
possible quantity that can show us an observable effect 
of chirality is the lattice current density. This operator is 
defined between two lattice sites indexed by two pairs of 
integers m and n, and situated at r m and r n . It measures 
the number of particles that flow from one site to another 
per unit time and takes the standard form [87] 



Jr 



i(^mn H~ ^nm)^n (^mn ^rim ) ^ in ■ (30) 



(a) 



(b) 




FIG. 5. (Color online) Local density of states: states are 
localized in the hexagonal plaquett es wh en the energy over- 
laps with the flat band at (f> = ~ (5(a) I; for the same sys- 
tem, intra-gap states are confined at the edges of the sample 
(5(b) l; dispersive middle band for (f> — J contains extended 



states (5(c) I; for comparison, local density of states for the 
disordered <f) = ^ system is shown in 5(d) for an energy E 
within the proximity of the middle bana! 



FIG. 6. (Color online) Plots of the expectation value of the 
lattice current operator in two implementations of the <f> = t 
system. Left panel: the system has a dispersive middle band 
overlapping in energy with the edge mode (See Fig. Ijjl. In the 
presence of a 5-function impurity the edge state current (red) 
will deviate around the impurity, maintaining its chirality, if 
the energy of the edge state lies in the gap and does not over- 
lap with the middle band. This corresponds to the anomalous 
quantum Hall phase. The current for a state whose energy 
is within the overlap region with the middle band is shown 
in green. Right panel: in the presence of a dispersive middle 
band (with (f> — ?) and disorder, the edge state can leak into 
the bulk with a finite probability. 



probability to scatter into the bulk, as shown in the right 
panelof Fig. [6j 



Using the notation of Sec. IV B the observable quan- 
tity of interest is the current expectation value for an 
eigenstate |^) of the tight-binding Hamiltonian J%? is 

<*|jm„|*>. 

From the study of edge mode dispersion relations in 
Sec. we have concluded that the introduction of a 
small impurity at the boundary can mix all momenta 
and the edge state can scatter into the bulk. We test this 
by introducing a <5-function impurity at one of the sites 
on the boundary of the system. The expectation value 
(^1 Jmn \&) is the familiar current density associated with 
the wavefunction 

We consider a system containing one impurity local- 
ized in real space at a single site. This impurity couples 
all pairs of k-points, implying that whenever the edge 
mode and the bulk state are at the same energy, but sep- 
arated in momentum, the edge mode will have a finite 
lifetime towards scattering into a bulk state. Results for 
the current operator in the Kagome system are plotted 
in Fig. [6] In the left panel, we have plotted the cur- 
rent of an edge state whose energy is located in the gap 
just above and not overlapping with the middle band at 
4> = j (see Fig. [I]) - the current (in red) will go around 
the ^-function impurity without leaking into the bulk. 
This corresponds to the anomalous quantum Hall phase 
and is reminiscent of the situation at <f> = As soon 
as the edge mode and the bulk band overlap in energy, 
current at the edge has a probability to scatter into the 
bulk. A state at this energy will essentially live inside of 
the bulk (current depicted in green). This corresponds 
to the anomalous Hall phase. If instead of an impurity 
at the edge we had strong phase disorder W$ = ^ for a 
system with flux ef> = ^, there would again be a finite 



BERRY'S PHASE AND ANOMALOUS HALL 
EFFECT OF NEUTRAL PARTICLES 



When a system in an anomalous Hall phase is placed in 
an external gauge field, it is possible to measure Berry's 
phases around closed loops in the Brillouin zone and band 
Chern numbers. This section is dedicated to such observ- 
ables. The treatment is applicable whenever the particles 
are neutral, and therefore suitable for cold atom systems 
as well. In Sec. |V A| we explain a scheme to realize syn- 
thetic gauge fields on the Kagome lattice that involves a 
time-modulated potential gradient in the cavity array, as 
suggested in Ref. [SHI In Sec. |VB| we present a way to 
measure the Berry phases of wavepackets of photons at 
a given energy based on interference. We also present a 
manner in which the Chern number can be directly ex- 
tracted from level-counting in the spectrum of the system 
in a uniform synthetic magnetic field in Sec. |V C| 



A. Realization of a synthetic gauge field 

We discuss here how to simulate the effect of a mag- 
netic field in a system of neutral particles. Such artificial 
magnetic fields, including frustrated or staggered config- 
urations, have been actively studied in recent years. In a 
cavity, photons coupled to artificial atoms, such as in the 
Rabi model, are subject to an effective gauge potential 
induced by the dipolar coupling [88] . On a lattice, gauge 
fields can in principle be simulated by creating effective 
Peierls phases through time-periodic driving as proposed 
in Refs. [54, 89, 90 for cold atoms. In what follows, we 
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shall follow the line of thought of Ref. [55] ■ The method 
relies on adjusting the potential terms individually for 
each site in the trapping optical lattice of the cold atom 
system. In cQED systems, this is equivalent to adjusting 
the frequency of each resonator individually and time- 
dependently, which has become experimentally possible 
[3H1 EH! ESI inU - R-ef- [35] has already proposed a means of 
realizing a synthetic gauge field in a square lattice array 
of coupled photonic cavities. 

Generally, a static gradient of the frequency w m of 
the resonator located at site m yields an artificial elec- 
tric field. Additionally, it is possible to have time- 
varying frequencies w mjT , which produce complex hop- 
ping amplitudes and so mimic the effect of a gauge field. 
Consider that that there is a driving frequency at site 
m = (mi, m-i) with the following dependence on the time 



<7G T — i/C 



m 



(T\uj + tvjjQ cos (f2r + m2#) mi) aJn a r 



(31) 

where the time-dependent Hamiltonian is the original 
tight-binding Hamiltonian of Eq. (J2j|, to which we 
add a time-dependent perturbation. The perturbation 
amounts to a tilt along the Ai direction, modulated in 
time. VL is the driving frequency, and for the following 
derivation we assume that the driving frequency is reso- 
nant with the on-site energy hio in the sense that f2 = uj. 
The details of this choice and the calculation are pre- 
sented in Appendix | A 1| Here we summarize the results. 

In the rotating wave approximation, the time depen- 
dence of Eq. (31) has the following effects: photons ac- 



quire an additional phase m-2,9 along each bond in the 
Ai direction; they acquire no additional phase along the 
bonds in the A2 direction. The following equations sum- 
marize the changes incured by the complex hopping am- 
plitudes: 



bond || Ai : |t|e _i * -> |t|e" 
bond || A 2 : |t|e _i * -> \t\e 



(32) 



The function J_\{z) is a Bessel function of the first kind. 
There will be a total phase f = 29 around the parallel- 
ogram plaquettes of area |Ai x A 2 | (see Fig. [I]). The 
artificial gauge field extracted from Eq. (31 1 gives rise to 



a uniform magnetic flux across the lattice, although the 
field itself is not uniform at the level of one unit cell due 
to the presence of the oblique bonds; nevertheless, this 
non-uniformity can be removed by a gauge transforma- 
tion (see Appendix All. This construction is equivalent 
to the Peierls substitution: a particle in a gauge field 
s/ s (y) acquires a phase 9 mn = ^ f* n dr.g/ s (r) between 
two sites located at r m and r n (the choice of path does 
not matter if the vector potential is that of an infinitely 
thin solenoid at the center of the plaquette). We shall 
make the convention that the photon coupling to such 
artificial gauge fields has an effective synthetic charge e s , 
and all fields to be used below will have a subscript "s" 



to highlight the distinction from actual electromagnetic 
fields. 

The experimental realization of an artificial field would 
allow us to measure Berry's phases around closed loops 
in the Brillouin zone or directly access the Chern num- 
ber, as explained in Sec. |V B| and |V C[ respectively. The 
potential to create a magnetic field can make a measure- 
ment of Landau levels of the Kagome lattice accessible, 
see Sec. IVCl 



B. Berry's phases from semiclassical dynamics 

In this section we describe a method to measure 
Berry's phases around closed constant energy contours 
in the Brillouin zone. For this we draw upon the semi- 
classical dynamics of a wavepacket within a given band 
of a Bloch Hamiltonian. 
In Sec. 



TVA" 



we have noted that, due to the po- 
larization phenomenon, wavepackets accelerated in the 
Aj direction drift in the A 2 direction. The motion of 
wavepackets subject to a uniform force is complicated due 
to Bloch oscillations. Proposals exist to access Berry's 
phases in a cold atom system relying on a measurement 
of the group velocity of a wavepacket under an external 
scalar potential exhibiting Bloch oscillations [53]. Here 
we present an alternative method to map Berry's phases 
of wavepackets, based solely on interference, and without 
the need to measure the group velocity of wavepackets. 

The semiclassical dynamics of a wavepacket centered 
at momentum and coordinate k c , r c and in a Bloch band 
is given by [92H94] 



k c x J? kc 



Kk c = -e s S s {r c ) - e s r c x 3$ s (r c ), (33) 

where $ s ,3@ s are synthetic classical fields, and e s is a syn- 
thetic effective charge, which is left as an explicit quantity 
in the equations in order to highlight the analogy with 
electronic systems. E(k c ) is the Bloch band energy, and 



= d k x #(k) 



(34) 



is the Berry curvature [Uj associated with the Berry 
gauge field introduced in Eq. ([7]); we have dropped the 
band index n, as we neglect transitions to other Bloch 
bands; & only has a z component. The semiclassical Eqs. 
( 33 1 are dual in the following sense: the second equation 



represents the electromagnetic force; in the first equation, 
the band energy plays the role of a scalar potential in 
momentum space, while the second term is nonvanishing 
only when the wavepacket is accelerated, and the Berry 
curvature acts like a magnetic field in momentum space. 
The group velocity is corrected by an "anomalous" 
contribution |95j which is now known to be inherently 
connected to geometric phases of wavefunctions j92j [96] . 

The geometric phase of a wavepacket can be measured 
in a system subject to a uniform synthetic magnetic field 
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perpendicular to the plane, 8§ s = \3§ s \t., but no elec- 
tric field. The semiclassical equations of motion for the 
wavepacket become 



— — k r X Z 



hk c = -^z, Js (k c )^- x <% s 

n ak r 



and Zgg, is a correction from the Berry curvature: 



Z&s ( k c) = 



1 



1 + 



(35) 
(36) 

(37) 



Eq. ( 36 1 implies that the force k c is always perpendicular 

dE 



to the group velocity gg-. Consequently, the energy of 
the wavepacket is a constant of motion, 



SE(k c 



dE 

dkr 



5k r = 0. 



(38) 



The trajectory in momentum space follows a contour of 
constant energy in the Brillouin zone. The trajectory 
in real-space described by Eq. (35), on the other hand, 



is related to the trajectory in momentum space by a % 
rotation and a rescaling. A wavepacket evolving on a 
constant energy contour will return to the original point 
in k-space. If the orbit is closed in the first Brillouin zone 
(continuous lines in Fig. [7]), then the trajectory in real 
space will be a closed curve as well. Alternatively, the 
orbit may be periodic in the extended Brillouin zone and 
follow a separatrix (dashed lines in Fig. [7|) , in which case 
the trajectory in real space will not be closed. Finally, 
the Berry curvature correction in Eq. ( 37 ) can only affect 
the rate at which the momentum k c varies, but not the 
form of the trajectory. The acceleration will be small in 
areas of strong Berry curvature, typically close to the K 
points. 

It is possible to determine resonance conditions for 
those paths which close in the Brillouin zone. Let us as- 
sume that we are in the case in which the wavepacket tra- 
verses a closed curve in phase space. Along a phase space 
loop of constant energy E, denoted C(E), the wavepacket 
will acquire a phase [92] [94] 



7c 



j> dk ■ M{k) + dr. (k - 



.si. 



(39) 



where the gauge field s/ s generates the magnetic field 
88 s = d T X sf s . This can be further simplified to 



7c 



r — — (r x SS a ) ■ r 



I dk ■ M{k) - [ dt 
Jc Jo 

Idk-M(k) + ~[ rft((rxJ,)T), (40) 
Jc 2h Jo 



where T is the period of the motion. The first term of 
the integrand is the Berry phase accumulated along the 




FIG. 7. (Color online) Constant energy trajectories in mo- 
mentum space, overlayed on the constant energy contours of 
the lowest band at cj> — ?. Trajectories can be closed within 
the first Brillouin zone resulting in closed real space trajec- 
tories or, if the initial momentum lies on the separatrix lines 
(see, for example, the dashed lines), then the trajectory in mo- 
mentum space can only close in the extended Brillouin zone, 
and in real space the particle traces a straight line. 



path, Tc = § c dkM(k). The second term can be written 
entirely in momentum space, 



■in 







It 

h 



dt( (f x 8§ s ) ■ r) 

h 



dt 



k x z x f , 



-k x 



2e s 



<j) dk'(kxz) 



(41) 



where dk • (k x z) = 2 dS? is twice the surface area ele- 
ment for the surface enclosed by the curve C(E) in mo- 
mentum space. Collecting the two pieces, the full phase 
acquired by the wavepacket around a closed loop in mo- 
mentum space is 



7c 



(42) 



The arriving wavepacket will interfere with the emitted 
wavepacket and the amplitude will be maximal if 



S^c = 2tt {n + u M ), 



integer, (43) 



where the Maslov index [97] is taken to be vm = h. 
There is an implicit energy dependence in the trajectory 
C = C(E), therefore the phases in Eq. (43) are energy 
dependent. This last equation is the Onsager relation 
[98j and can be viewed as a Bohr-Sommcrfcld quantiza- 
tion condition for the closed trajectories in phase space. 
An analogous situation occurs in electronic systems, in 
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the de Haas - van Alphen effect [§3], where peaks in the 
magnetization as a function of the external magnetic field 
are a result of such a resonant behavior. 

In a photon system, the Berry phases can be accessed 
as follows. Suppose that a wavepacket of known energy 
E and initial momentum k c (i = 0) can be produced. 
Upon tracing a closed path in real space, the traveling 
wavepacket returns and interferes with the emitted sig- 
nal. One can tune the external field S8 S such that the 
amplitude at the emitter is resonant. Two consecutive 
resonances determine the area S?c of the common closed 
path in momentum space, see Eq. ( 43 1 . This allows FIG 



for the determination of the Berry phase along the path 
C(E) up to a multiple of 2tt. By keeping the energy E 
fixed and changing initial momentum one can explore all 
curves at the given energy within the Brioullin zone. Two 
examples are given in Fig. [7J What such a measurement 
would yield with and without disorder is presented in 
Sec. I VII in connection to a realization of the anomalous 
Hall effect in the Kagome photon system. 

We have shown that an interference experiment can be 
realized to measure the Berry phase of a wavepacket in- 
jected at a specific energy and initial momentum. Berry 
phases appear in measurements of a photonic equivalent 
of the anomalous Hall effect, as explained in the next 
section. 

Moreover, the measurement can be changed to probe 
Landau levels. Instead of keeping the energy fixed and 
varying the magnetic field 3& ' a> let us keep the field fixed 
and vary the energy. Separations between resonant or- 
bits discussed above correspond to transitions between 
Landau levels. The analysis of level splitting in the pres- 
ence of 38 a allows us to directly probe the Chern number 
of a Bloch band. This is the subject of Sec. |V C| 



C. Direct determination of Chern number from 
level counting in the Hofstadter spectrum 

In this section we present a way to directly determine 
the Chern number of a Bloch band from the spectrum in 
a synthetic magnetic field. In a magnetic field, the three 
original Bloch bands will split into subbands, sometimes 
called magnetic Bloch bands. The resonant semiclassi- 
cal trajectories of Sec. |VB| correspond to the magnetic 
subbands of the original Bloch band under the influence 
of the synthetic magnetic field 38 a . The Chern number 
influences how a Bloch band splits into magnetic sub- 
bands 92 . There exists a maximum number of closed 
resonant trajectories within the Brillouin zone. This is 
obtained by observing that the maximal area of a single 
resonant closed trajectory C, has to be equal to the 
area of the first Brillouin zone. This dictates, then, that 
the number of subbands obtained by splitting any one of 
the Bloch bands is given by: 




Low magnetic field SS S spectrum of the Kagome 
lattice with (f> = j : at external field, we recover the original 
bandstructure of Fig. [I] At flux p/q there are 3q subbands 
with possible degeneracies. For example, at p/q = 1/5, the 
density of states reveals that there are 6 separated subbands in 
the upper band, 5 degenerate subbands in the middle band, 
and 4 separated subbands in the lower band in agreement 
with Eq. ( |44[ ), and consistent with the fact that the lower 
and upper Bloch bands have Chern numbers -1 and 1, respec- 
tively. The bands become degenerate (no energy gaps in the 
spectrum) at p/q = 1/4. See text for details. 



where v is the Chern number associated with the band, 
/ = ,g?3 ° ^ is the flux through the unit cell of a uniform 
synthetic magnetic field 38 s , in units of ^, and the square 
brackets indicate the integer part of the real number. The 
field 38 s has to be small enough such that the resulting 
subbands cannot become degenerate. This is argued in 
more detail numerically in what follows. 

We are aiming to solve the Hofstadter problem of a 
tight-binding system in a magnetic field [73] . Let us pick 
the Landau gauge, srf a — (—3§ s y, 0), where the two com- 
ponents are cartesian. Due to the presence of a magnetic 
field, the Hamiltonian will couple different points in the 
Brillouin zone. If the flux / introduced above for the uni- 
form synthetic field is a rational number expressed as p/q 
with p and q non-negative integers which are relatively 
prime, then, in general, the resulting Hamiltonian will 
couple k y to k y ± ^£ and to k y ± . This momentum 
space coupling can be removed by remarking that trans- 
lational invariance is recovered if one reverts to a q times 
larger unit cell in real space, and a q times smaller Bril- 
louin zone, called the magnetic Brillouin zone, defined 



for our lattice as 



0. 



qn 



-^ ZL - 



The original 3-band 



D = 



(44) 



Hamiltonian becomes a 3q-band problem defined on the 
reduced Brillouin zone. Folding of bands appears, and at 
a general level the 3q bands are allowed to become de- 
generate within the magnetic Brillouin zone. Technical 
details are given in Appendix |A 2| For now we focus on 
counting the magnetic subbands obtained from a given 
Bloch band, which provides a direct way of determining 
the Chern number of the original (zero field) bands. 

Assume for simplicity that the dimensionless flux per 
unit cell is / = g for some positive integer Q. For large 
enough Q, the field is weak, and any one of the three 
original Bloch bands will split, according to Eq. (44 1, 
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into [v + Q] , where v is the Chern number on each of the 
three original Bloch bands. This is consistent with the 
fact that the spectrum in a magnetic field at flux 1/Q 
should exhibit 3Q levels. Taking the band-structure at 
4> = j (see Fig. [Tj) as representative for our time-reversal 
symmetry broken phase with gapless edge modes, the 
formula in Eq. (44 1 predicts that the lower, middle and 
upper bands will split into Q — 1, Q and Q + l subbands, 
respectively. The caveat to this discussion is that for 
special fractions / one can expect that the resulting sub- 
bands become degenerate even if the original system is 
non-degenerate [T5"l l9"2"] ■ 

One can directly access the Chern number of a spe- 
cific band as follows (see Fig. [8]): set the external syn- 
thetic field to a specific rational number dimensionless 
flux and count the resonant trajectories corresponding 
to each of the three original Bloch bands (they will be 
typically separated by gaps of order \t\, as in the original 
system, if the field is weak) . For example (see Fig. [8]) , at 
/ = 1/5, the lower, middle, and upper Bloch bands split 
into [v+5] = 4, 5, 6 subbands, consistent with their Chern 
numbers v = —1,0,1, respectively. There are a total of 
15 subbands, but the 5 subbands corresponding to the 
middle band are degenerate. Importantly, at / = 1/4, 
the spectrum has no energy gaps, in which case Chern 
numbers are exchanged between bands. We have proved 
numerically the requirement / < 1/4 to probe the Chern 
numbers of the original Bloch bands. Since the spectrum 
of the Hofstadter problem is periodic in /, this require- 
ment actually translates to / € [0, 1/4) U [2, 9/4) U .... 



D. Bounds on the strength of the synthetic 
magnetic field 

When the external field is too strong, or varies too 
quickly, tunneling to a different band becomes possible. 
To be safe from Landau-Zener tunneling [OS], the mag- 
netic field has to be small. This translates to a condition 
[94] on the period of motion on a closed loop, T, intro- 
duced in our calculation above, 



n „ Ie, 



(45) 



where E g is the size of the gap, which in our problem is of 
order the hopping energy, E g ~ \t\, and E is the constant 
energy along the loop C(E). The period of motion can 
be further reexpressed by making an estimate from Eq. 

PL 



T 



hlc(E) 



(46) 



where the averages are taken over the path C, whose 
length in momentum space is denoted lc(E). Collecting 



equations yields the following condition 



hl c {E) 




(47) 



It turns out that this bound is independent of the lattice 
spacing a or the absolute value of the hopping integral 
\t\. For the typical values of the Berry curvature etc., the 
bound of Eq. ( 47 ) is satisfied by taking the flux of £3 S , 



f « 1, which is consistent with the results in the Sec. 

EH 



VI. ANOMALOUS HALL EFFECTS OF 
PHOTONS 

In this section, we discuss the anomalous Hall effects 
of light and evaluations of Chern numbers for the clean 
and for the disordered case. 



Chern number 



The experiment of Sec. |VB| amounts to measuring line 
integrals of the Berry gauge field ,S?(k) [IB] : 



T C (E) 



j> rfk.^(k). 



(48) 



Since C(E) is a constant energy curve in the Brillouin 
zone, the phase can be reexpressed as a sum only over 
the states of the Bloch band which lie below the given 
energy of the curve, E, 



d 2 kd(E - £(k))jF k . 



(49) 



BZ 



This value is, in general, not quantized: only after sum- 
ming over the full Bloch band can the phase become a 
multiple integer of 2n. In electronic systems, the non- 
quantized part is the intrinsic contribution to the anoma- 
lous Hall effect [Ml 1100] of a partially filled band (non- 
quantized Hall conductivity cr xy ), and can be interpreted 
as a Berry phase of quasiparticles at the Fermi surface 
[51] . In a bosonic system, we interpret this as the anoma- 
lous transport of a wavepacket whose energy overlaps 
with that of a bulk band. 

So far, we have dealt with a single Bloch band. In a 
multiple band system, we can define the following quan- 
tity as a sum over states below some energy: 



m J BZ 



(50) 



where now the Berry curvature is defined as in Eqs 



(J7J) and (|34|) but for the n th Bloch band \nk). v{E) is 
plotted in Fig. [9j When the energy E lies in a gap, v{E) 
is an integer, and it is the sum of the Chern numbers of 
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FIG. 9. (Color online) Panel A: Dispersive middle band for tf> = ~. The Chern number, computed with formula of Eq. (501 
gives a non-quantized value in the overlap region which can be measured directly in an interference experiment. Panel B: Chern 
number calculation for the flat band system <f> = ^ with increasing amplitude of phase disorder: W$ = 0, f'fifif'f - Panel 
C: flat band system with increasing amplitude of on-site disorder W S it e = 0, J=L, \t\, The disorder potential has no spatial 
correlation, and is drawn from a uniform distribution (see text). 



all bands that lie below the given energy. We may rewrite 
this quantity in a form that is manifestly invariant under 
gauge transformations of Bloch vectors, and will be useful 
as we tackle the disordered case in Sec. IVIBI We define 
the projector P^ = Pk(E) onto the Bloch states of energy 
below E. The projection operator is obtained from the 
Bloch eigenvectors of the tight-binding problem, and is 
therefore k-dependent. We can rewrite Eq. (50 1 as 



v{E) = ^.Tr {P k [d kx P k , d ky P k ] } . (51) 

The number v{E) is an integer only when the energy lies 
between two bulk bands without touching them, inside 
of a gap, and it takes non-universal values for the slight- 
est overlap in energy with the bulk bands (see Fig. [9]). 
The quantization of this quantity corresponds to the ex- 
istence of counter-propagating edge modes (see Sec. Ill I. 
The detuning of the phase cf> will make v(E) take a non- 
quantized value, which corresponds to edge modes that 
can decay into the bulk of the material. We corroborate 
that disorder is similar in effect to changes of the phase </>: 
it also spreads the middle band leading to non-quantized 
values of the Chern number. This is explained in the 
next section and summarized in Fig. [9j 



B. Disordered system 

In this section we investigate the robustness of edge 
states in the disordered system. The Bloch state formu- 
lation used so far cannot be used since the system breaks 
translational invariance. To arrive at the formula for the 



Chern number in real-space for disordered systems, we 
start from the clean case. The trace in Eq. 



(51 



can be 

taken with respect to the basis of states localized at each 
site, |m), where m = (mi, m?) indexes a site on the lat- 
tice. In this section we return to the x, y basis, as this is 
more useful for deriving formulae suitable for numerical 
computations. Letting as usual r m = miAi + ra-i A2, 



v{E) 



lim 

N->oo 



2ni 



^(E) [-i[x,P(E)],-i[y,P(E)]] 



(52) 



where P(E) is the Fourier transform of the operator 
Pk{E) defined above, and projects onto states with en- 
ergy below E. The total number of sites in the system is 
denoted by N. In the thermodynamic limit for a clean 
system, this formula is exactly that from Eqs. (50) and 
(51). Its advantage is that it allows us to compute the 



Chern number in a disordered system. In Ref. |101j it 
was shown that the Chern number v(E) is an integer as 
long as the energy E belongs to a gap in the spectrum. 
The integer value can change to a different one only if 
the energy E crosses a region of extended states, such as 
traversing a bulk band. 

When the system is disordered and finite, the prob- 
lem of averaging over disorder configurations arises along 
with that of taking the thermodynamic limit. The Chern 
number can be defined as an ensemble average over dis- 
order configurations |101l 1102] 



(53) 



v{E) = / dfi(6)u s (E), 
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where the disorder configurations S are distributed ac- 
cording to a measure dfx(6), and v$(E) is the value of 
the Chern number for the specific disorder configuration. 
One defines a Hamiltonian J£§ for each disorder config- 
uration S, and associates the projector Ps(E) to it as in 
the clean case, which gives a formula for v$(E): 

ME) = 

-^T.^W) [-i[x,P6(E)]-i[y,P 5 (E)}}\r m ), 

m 

(54) 

where now x, y are the position operators in cartesian 
coordinates. This equation has the property that the 
average over an ensemble of disorder configurations can 
replace the thermodynamic limit of a single system. The 
technical details, and reformulations suitable for numer- 
ics are reserved for the Appendix |B) 

We have computed the Chern number for the Kagome 
system with a flat middle band using the methods de- 
scribed in [102] at (f> = |, implemented on a lattice of 
24 x 24 sites, or 192 unit cells. At each point, the Chern 
number was computed for 40 disorder configurations and 
averaged. A white noise disorder potential (uncorrelated 
from site to site) was sampled on the lattice. For a disor- 
der amplitude W, a random number between [—■?-, ^-1 
is produced. The two types of disorder that can appear 
on the lattice are scalar disorder, on the on-site frequen- 
cies, which we denote W S it e , and vector disorder on the 
hopping phases, denoted W^. 

Fig. [9]shows our results for W S i te between and 1.5|i| 
and Wrf, between and The originally flat middle 
band spreads with disorder; in view of Sec. |IVB[ the 
broadening of the middle band is associated with states 
leaking out of the hexagonal plaquettes due to the de- 
tuning of the flux. Comparatively, the role of scalar or 
vector disorder potentials is the same. We are not ad- 



dressing here the interesting question of comparing the 
two types of disorder. 

VII. ANISOTROPIES 

Recently it has become possible to realize, shift, 
and merge Dirac cones in optical lattices [1031 H04j . 
Anisotropics on the honeycomb lattice have been stud- 
ied [105] . In the cQED realization of the Kagome lattice, 
anisotropies are inherent: the on-site energies Huja.b,c 
and hopping amplitudes t A B.c obtained after integrat- 
ing passive circuit elements are generally distinct (this is 
presented in detail in Ref. [21]). Small anisotropies leave 
the main features of the bandstructure intact. Below is 
a quantiative analysis of the effects of anisotropies. 

Dirac cones are stable due to discrete symmetries 
[1061 1107] . The isotropic Kagome lattice (isotropic hop- 
ping amplitudes and site potentials), is symmetric under 
inversions with respect to the center of a hexagonal pla- 
qucttc, and under rotations by 27r/3, up to the permu- 
tation of the sublattice indices (lattice geometry in Fig. 
[lj. Additionally, for special values of the phase <f> the 
system is invariant under time-reversal symmetry. These 
three symmetries fix the position of the Dirac cones at 
the corners of the Brillouin zone between two bands. Two 
more Dirac cones are merged at the T point (where the 
dispersion of one of the bands is flat, and of the other 
quadratic). If the rotational symmetry is broken due 
to anisotropy, then the Dirac cones can be shifted in k- 
space. For strong enough anisotropy, Dirac cones meet 
at time-reversal invariant points in the Brillouin zone, 
where they can annihilate. Such points are the T point 
and the M points. In the example below, Dirac points 
annihilate at the M points. 

Using the notations introduced in Sec. [TTJ let us de- 
fine 01,2,12 = cos(ai.2,i2)- The Hamiltonian including 
anisotropies is 



ftjJJA 2t A B COS Ct\ 2t A cCOSa 2 

= | 2tBACOSa\ Hujb 2±bc cos Ct 12 

2tcACOsa 2 2tcscosai2 fiwc 



(55) 



which leads to the following equation for the three energy levels: 

- E 3 + h{oj A + lob+ lo c )E 2 - (h 2 (uj A ujB + ujbojc + u c u A ) - ^{\t AB \ 2 c\ + \t AC \ 2 c 2 2 + \t B c\ 2 c 2 12 )) E 

-Ah(uj A cl 2 \tBc\ 2 + u>Bc\\t AC \ 2 + io c c 2 \t AB \ 2 ) + l6\t AB t Ac t B c\cic 2 ci 2 cos(3</>) + h 3 u A u) B uc = 0. (56) 



Let us consider, for example, the case when the hopping 
integrals are isotropic, but there is an anisotropy between 
the on-site energies hoj Al B,c at the level of each triangular 
plaquette. The spectrum can become degenerate at the 
M points defined by k x = 0, k y = ±^tj- One has c\ = 1, 
c 2 = 0, and C\ 2 = 0. Then for all 4>, a degeneracy occurs 



r 



if = (u A + lob ± \J Ju~a - ^b) 2 + P]*ab| 
the other two pairs of M points, k x = ±?, k 



/2. At 



« 



± 



~^aV3 d e g enerac y occurs if the 
cyclic permutations of the condition above hold. This 
procedure allows us to produce a degeneracy between 
the upper and middle bands or the lower and middle 
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bands (although all three bands cannot be degerate at 
the same k). The values for uja,b.c necessary for this 
degeneracy are on the order of the hopping amplitude 
\t\. The touching cannot be lifted by varying the phase </>. 
The degeneracy at the M points is quadratic in at least 
one direction in k-space, which is a feature of the fact 
that two Dirac cones with linear dispersion annihilate at 
the time-reversal invariant point (see, for example Ref. 
106 ). One can tune the anisotropy as to annihilate a 
pair of Dirac cones between the lower and middle band, 
for example. As soon as the Dirac cones have annihilated, 
the bands are no longer degerate and the respective gap 
no longer carries edge modes. 

In an experiment it is plausible to consider the case 
of anisotropic phases 3</> np and 3</i 'down around the up- 
pointing and down-pointing triangular plaquettes, re- 
spectively. This arrangement maintains the requirement 
that there is zero flux per unit cell, and the effect will be 
to merely shift the Dirac cones due to the anisotropy. 



VIII. INTERACTIONS 

In this section we address the stability of the topolog- 
ical phase in the presence of interaction effects. A host 
of methods have been used to show that the topologi- 
cal phases are stable even in the presence of interactions 
[108| I109j . At a simple perturbative level, we show that 
interactions cannot change the effective free particle dis- 
persion bandstructure: they cannot cause the bands to 
become degenerate and reopen the gap. The general ar- 
gument goes as follows: if the system is in a topological 
phase with gaps between bulk bands, then it is enough to 
show that for weak enough interactions the shifts in the 
energy levels are smaller than the free particle spectrum 
gap size of order the hopping strength \t\. 

The following are interaction Hamiltonians which can 
be potentially implemented in arrays of resonators. The 
boson Hubbard model was discussed, for example, in an 
important paper by M. P. A. Fisher et al. |110j . The 
Jaynes-Cummings Hamiltonian [81] . on the other hand, 
models an interaction between lattice photons and quan- 
tum two-level systems situated at each site. Both the 
boson Hubbard model and the Jaynes-Cummings model 
exhibit a phase transition from a superfluid to a bosonic 
Mott insulator state. Both interactions are in principle 
realizable in a cQED experiment [24l [26] . For complete- 
ness, we add a brief analysis of fermionic interactions in 
Appendix [C| 



A. Bose-Hubbard model 

We start with the unperturbed Hamiltonian, the one 
presented in Eq. ([3]), 

* = E V^kV'k = £ aU-aZTVk, (57) 

keBZ a,^,k<EBZ 



where the Greek indices indicate the sublattice in the 



spinor = (a^. 



Ak> a Bk> "Ck 



The interaction Hamilto- 



nian for the Bose-Hubbard model is quartic in the cre- 
ation and annihilation operators 



J^bh = -7T E o) m a m {a) m a m 



I) 




^ a a,k 1 -k a Q,k 2 +k a a>k 2 a Q!,ki , 



k!,k 2 



(58) 



where TV is the number of unit cells and we have used 
the Fourier transform on each sublattice a = A, B, C as 
a Q k = J2 m e a e~ lkR "" a m- We would like to describe 
the spectrum of elementary excitations above the ground 
state. The minimum of the spectrum is at the T point: 
the ground state of the bosonic system has all photons 
condensed at k = in the lowest band (see Fig. [I]). Let 
us denote the minimum single particle energy at the T 
point by Eq. 

In the following treatment we shall restrict to a sub- 
space of constant particle number n, where n = Nuq. 
We have denoted by N the number of unit cells in the 
lattice, implying that there are 37V sites, and by n the 
number of condensate particles per unit cell. The spec- 
trum of elementary excitations in the superfluid phase is 
determined by using the Bogoliubov approximation |lllj , 



yj Nn a 8 v 



(59) 



Eq. ( [59] ) defines the excited state operators b a ^ {) . The 
full Hamiltonian can be brought into the following form 



ka/3 



Un 
2 



E 



fra,k&a.-k + b\ v b] 



a,k a, — k J ' 



k/0,a 



,(60) 



where Eq = nE is the total ground state energy, and 
the remainder describes the excitation spectrum. We 
shall only focus on the excitation spectrum and drop the 
ground state energy from our notation. The Hamiltonian 
can be diagonalized by a Bogoliubov transformation [lllj 
from the particle operators b a ^ to a new set of quasipar- 
ticle operators &„ ; k, which annihilate a quasiparticle in 
the n th band. The details of the Bogoliubov transforma- 
tion are reserved for Appendix [D] The Hamiltonian of 
Eq. ( 60 1 can be recast into a diagonal form from which 



one can extract the quasiparticle dispersions, which we 
denote by £ n (k), 



k,n 



(61) 
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Here (7^ = cTx.m±itTy,m are Pauli ladder operators 



FIG. 10. (Color online) The left panel shows the deviation of 
the effective quasiparticle dispersion f«(k) (black) from the 
free particle dispersion E n (k) (dashed, red) for the <j> — ^ 
system with U/\t\ — 0.13 (the plot is along the line k x = 
0). Right panel: low energy quasiparticle spectrum exhibits 
linear dispersion (magenta) at the T point k = 0; plotted for 
comparison is the free particle dispersion (light gray). 



In the following, we discuss the stability of the topo- 
logical phase when the Bosc-Hubbard interaction is in- 
troduced. In general, if the interactions induce no band 
crossings, the Chern numbers on individual bands will 
be conserved, and the topological phase is maintained. 
The Hubbard interaction induced dispersion £„(k) can 
only change significantly from the free particle dispersion 
E n (k) of Eq. ([5| in the vicinity of the band minimum at 
the T point, where the quasiparticle has a linear disper- 
sion (the sound mode). Fig. 10 shows the quasiparticle 



dispersion in comparison to the free particle dispersion, 
in the flat band system at <f> = ^. As long as the Bose- 
Hubbard coupling U is on the order of \t\ or smaller, the 
shift in the energy levels cannot make the gap close and 
reopen, and the topological phase is protected. In the 
weak interaction regime, the middle flat band can ac- 
quire a dispersion: This situation is similar to the effect 
of detuning the hopping phase <f>, or the effect of disorder. 



B. Jaynes-Cummings interaction 

We proceed with the Jaynes-Cummings interaction [ST] 
which couples photons to quantum two-level systems 
at each site. The Jaynes-Cummings system exhibits a 
quantum phase transition of polaritons from a superfluid 
phase to a Mott insulating phase [35] , resembling that of 
the Bose-Hubbard model. The system is governed by the 
following Hamiltonian which describes the tight-binding 
photons, the two-level systems, and the coupling between 
these 

= X! *m,n a L a n + 2J eO"mCr m + ^ 9(&m a m + h.C.) 
m,n m m 

(62) 

where we have written the tight-binding Hamiltonian in a 
more generic form as a sum over all pairs of sites, in terms 
of generic hopping integrals t m n . Note that here the on- 
site energies Huj have been included in the diagonal terms 



describing a two-level system at lattice coordinate m. 

To characterize the stability of the anomalous Hall 
phases in the presence of a weak (finite) Jaynes- 
Cummings interaction, in Appendix [E] we find an effec- 
tive low-energy Hamiltonian for the photon sector Jf^. 
We consider the following limit of small conversion g, 



g < |t| < fvjJ. 



(63) 



This assumption is prompted by the following stability 
condition specific to the Jaynes-Cummings lattice [33] . If 
the ground state energy of the single particle spectrum, 
Eq, has a negative value, there is no mechanism to limit 
the number of photons that populate the ground state, 
and therefore the system becomes unstable. We there- 
fore need Eq > 0, which translates to a condition that 
Hut > z\t\, where here z = 4 is the coordination number 
of the Kagome lattice. In a typical cQED experiment, 
we expect \t\ to not exceed 100 MHz and lu to lie in the 
GHz range. Additionally, under reasonable experimen- 
tal conditions the two level system excitation energy e 
and the (renormalized) resonator frequencies u) will both 
lie in the (microwave) GHz range. The condition that 
Hu> « e is sufficient to order the two-level systems in 
the ground state, and the determination of an effective 
photon-photon interaction is well motivated. This, as 
we show below, leads to a positive value of the Hubbard 
coupling, i.e. a repulsive interaction. 

To fourth order in the Jaynes-Cummings coupling g, 
wo obtain the following effective low-energy Hamiltonian: 



JT eff = 

ha 1 

EQ-t 



E {] 



a lo a w 



f49 4 
Wo - t 



-^a w ai a lo ai , 



(64) 



where fa and fa are positive dimensionless factors which 
can be found in Appendix [e| and the operator a\ Q creates 
a particle in the ground state of the photon tight-binding 
Hamiltonian (lower band at the T point, k = 0). Eq. 
( 64 1 shows that to lowest order the interaction shifts the 



(lowest) band minimum by an amount j^ 9 _ t ■ The next 
non-vanishing contribution appears at tne quartic or- 
der. The repulsive character is obtainable under reason- 
able experimental conditions for the Jaynes-Cummings 
model, where both Hlo and e would be in the microwave 
range. One way to understand the presence of a repul- 
sive photon-photon interaction is by noting that the two- 
level systems are bosons with hard-core on-site repul- 
sion. More precisely, if two-level systems are represented 
by bosons with interactions, the Jaynes-Cummings term 
can be absorbed by a shift of the bosonic two-level sys- 
tem field, and then the interaction term of the two-level 
systems also generates an on-site interaction for the pho- 
ton in the shifted basis. The correspondence between the 
Jaynes-Cummings lattice and Bose-Hubbard model has 
also been demonstrated through a field theory approach 
close to the Mott-superfluid transition of light [33 . In the 
calculation that led to Eq. (64), we have considered the 



17 



ground state with all photons at wavevector k = 0. This 
does not provide precise information about the range of 
the interaction, but it is telling that a Hubbard-type re- 
pulsive photon-photon interaction can emerge. Finally, 
as long as the perturbation g is weak compared to the 
size of the gap, which is of order \t\, the detuning e, and 
the on-site energy hui, the interactions cannot cause de- 
generacies, and the topological phase will be stable. 



Appendix A: Synthetic gauge fields 

In this appendix we present a more detailed derivation 
of the results in Sec. IV Al 

1. Magnetic field 



IX. CONCLUSIONS 

Motivated by the recent experimental progress in the 
context of arrays of electromagnetic superconducting res- 
onators jinj ! we have investigated the anomalous Hall ef- 
fect of light on the Kagome lattice with artificial gauge 
fields [21] . The photonic system here exhibits equivalents 
of the quantum Hall effect without Landau levels, and 
the anomalous Hall effect with a non-quantized Chern 
number. In particular, we have shown that a topologi- 
cally trivial band can affect the quantization of Chern 
numbers as well as the robustness of the chiral edge 
modes. We have discussed observables which are acces- 
sible experimentally. We have introduced a method to 
measure Berry's phases around loops of constant energy 
in the Brillouin zone. The method is based solely on 
wavepacket interference and can be used to determine 
band Chern numbers or the photonic equivalent of the 
anomalous Hall response. It provides an alternative to a 
recent method proposed to measure line integrals of the 
Berry gauge field in cold-atomic systems, which relies 
on the measurement of group velocities of wavepackets 
and a force-reversal protocol [53 . In cQED systems, this 
is realized by adjusting the frequency of each resonator 
individually, which has become experimentally possible 
[3"5l IS"5l 151)1 |9"T] . It shall be noted that interference exper- 
iments can also be envisioned to probe the Landau levels, 
emerging when placing the Kagome lattice in a uniform 
magnetic field. An open and interesting question to in- 
vestigate in the near future would be the influence of ar- 
tificial gauge fields on the superfluid-Mott transition of 
light in cQED photon-based lattices, following a similar 
line of thought to that used in recent theoretical investi- 
gations |109L I112j . Recent progress in this direction has 
been realized in Ref. j!13j where a chiral Mott insulator 
with a gap to all excitations and staggered fluxes has been 
found. Transport of microwave photons through cQED 
lattices with artificial gauge fields and disorder would be 
an interesting topic to explore. Another relevant subject 
to consider more thoroughly both theoretically and ex- 
perimentally is the realization of quantum impurity mod- 
els exhibiting fingerprints of many-body physics such as 
the Kondo model and Spin-Boson models [114H117] . 

We acknowledge useful discussions with Stephan 
Rachel, Peter Orth and Arun Paramekanti. AP and 
KLH are supported by the NSF through DMR-0803200 
and DMR-0653377. AAH acknowledges support from the 
NSF Grant DMR-0953475. 



This derivation starts from an idea used in Ref. |89j 
for square lattices. Consider the time-dependent Hamil- 
tonian of Eq. (31 ), which we reproduce here: 



Jtf? T = iff + (hoj + hu>o cos (Qt + m20) mi) a>in a m- 

m 

(Al) 

We construct the solution to the Schrodinger equation 
as follows. If the hopping is supressed, \t\ = 0, then the 
following function solves the time-dependent Schrodinger 
equation: 



w 



<Pm{ T ) — ~ UJT — q sin (Or + Bm-z) rn\. (A2) 
The solution for \t\ ^ is constructed from this as 

|V)=5> m e^"|m), (A3) 



where the d m must now obey the following differential 
equation 



ihd r 



/ ' m.ri ' «*n 



(A4) 



In this last equation, the i mn is the tight-binding hop- 
ping integral, which takes the value |i|e ±l< ^ for nearest- 
neighbors and zero otherwise, and 4>m are the phases 



computed in Eq. ( A2 1 . The perturbation has induced 
a time-dependent phase factor e l ^ m ~^^ which we will 
now simplify by keeping only those parts that oscillate 
very slowly (rotating wave approximation) . We shall use 
the following expansion 



%z sm a 



£ 

I — — oo 



e lla Ji(z), 



(A5) 



where Jiiz) are the Bessel functions of the first kind. 
Upon inspection of the expansion in Eq. ( A5 ) and of the 
phase 4> m in Eq. (A2) we find that, in general, one must 



have the driving frequency Q be an integer multiple of the 
on-site frequency w, i.e. Ilo, in order to obtain at least 
one time-independent term in the expansion. The largest 
contribution is obtained if we take — uj. We obtain the 
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following effective changes to the hopping amplitudes 



bond || Ai : 

\t\e~** -> |t|e -< * 

bond || A 2 : 

\t\e-^ -> \t\e-^ 



bond || 

|t|e"* 



(A 2 -Ai) 

\t\e-^ 



J-i 



J- 



K + 1)) e 



(A6) 

The ratio ^ provides an additional experimental param- 
eter to tune the hopping strength via the Bessel functions 
of the first kind J-\. 

The time-dependent perturbation has induced spa- 
tially dependent phases and a dressing of the hopping 
integral \t\. These spatially dependent phases mimic the 
phases that would be produced by a gauge field in the 
minimal substitution. The fact that an additional phase 
is acquired along oblique bonds parallel to A 2 — Ai im- 
plies that at the level of each unit cell the phases will cor- 
respond to a field that is non-uniform across the unit cell. 
However, the total phase acquired by a photon traversing 
around a parallelogram unit cell of area |Ax x A 2 |, see 
Fig. [T] is going to be a constant equal to f — 29. Since 
the uniformity of the field at the level of the unit cell can 
be recovered by a gauge transformation, we perform all 
of our calculations (Appendix A 2 ) for a uniform field. 



2. Spectrum in a magnetic field 

In this appendix we show the detailed calculations for 
the spectrum of the Kagome system placed in a uni- 



form magnetic field SS a . The phase acquired by a par- 
ticle along an elementary parallelogram unit cell of the 
Kagome lattice (see Fig. 1} is / = = 8tt/ = 

•\/36a, where / is dimensionless, and b has units of in- 
verse length. The phase acquired by a photon around the 
unit cell is eight times the phase acquired on a counter- 
clockwise loop around a triangular plaquette. Let us for 
simplicity pick the following gauge field in the Landau 
gauge, g/ s — (—& s y,0), where the two components are 
cartesian. Due to the presence of a magnetic field, the 
Hamiltonian will couple k y to k y ± | and to k y ± b. With- 
out loss of generality, we take Hu> = 0. The Hamiltonian 
reads 



' ^ e ' k ' Al ' a Ak x ,k y a Bk x ,k y +b + 0>Ak X ,ky a Bk X ,ky-b 



+e lf s e 



j_ 



,-ik.A 2 t 



C . k x .ky B ,k x ,ky ~\~ r 



a 



C,k x ,ky B jk x jky , 



,-ik.(Ai-A 2 )-i 



+e "is i e 

+ e 8 a A,k x ,K a C,k x ,k v 



baV3 f 

^4. jkrj- jky k^ yky ~\~ T 

f h.c. 



(A7) 



For rational values of / = | , where p, q are relatively 
prime integers, one can reduce the Brillouin zone from 



the original [0, 



0, 



2tt 



0. 



4-7T 



a\/3 



o, 



a\/3 



to the magnetic Brillouin 



The couplings between different 

momenta disappear and we have replaced the original 
problem with that of a periodic one-dimensional chain of 
3q sites. Let us take a generic wavefunction to be 



9-1 

\ip) = ^2 Vwa 

n=0,a=A,B,C 



t 



10), 



(A8) 



where a) 



, „ h creates a photon on sublattice a at 

a given momentum. Then the Schrodinger equation is 
equivalent to the following set of three Harper equations 



Ek x ,k°i>Am = |t|e^-*A (e+^^cm-i + e-< k *> k >™i)^-^i> c , m+ x) + |*|e-^(e-<^+^)- A ^B,m+2 + ^, m - 2 ) 
E kx , k oJ Bm = \t\e^(e< k ^ m - 2 ^)- A ^ A , m . 2 + VA, m+2 ) + \t\e-«-&(A i 'W m -Vl)- A >1,c, m -i +^, ro+ i) 



B,m-1 



+ |t| e -^+i^^ e <(fc«,fc5+("»-l)|)-(Ai-A 2 )+i6^^ 



_ i baV3 

A,m-1 + e 8 WA,m+l 



(A9) 



The solution to these equations gives the spectrum of the system for every rational flux / = p/q, and the pattern of 
splittings into magnetic subbands is known as the Hofstadter butterfly |73j . The spectrum of the problem is periodic 
in / of period 2, for example f — and / = 2 systems have the same spectrum etc. 

I 



Appendix B: Evaluation of Chern numbers in 
disordered systems 



of Eqs. ( |53| and ( 54 ). The general theory was introduced 
in |101l 1102) . To keep our discussion generic and more 



This appendix is dedicated to presenting more rigorous 
statements behind the real-space Chern number formulae 
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suitable for numerics, we shall rescale our Brillouin zone 
to [0, ?f] x [0, ^] . The results for the reciprocal unit cell 
of the Kagome lattice can be adapted from the following 
by the introduction of the appropriate Jacobian. 

Consider a two-dimensional lattice with K orbitals per 
site. We may take {|ma)|a = 1, ...,K} to be kets local- 
ized at each site m = (m 1 ,m 2 )j corresponding to each 
orbital. This amounts to a basis of the full Hilbert space. 
The orbitals a may represent actual orbitals, or spin, 
isospin due to more ions per unit cell etc. The Bloch 
transformation is a unitary transformation that takes the 
Hilbert space H to a direct sum of C K spaces of if-tuplcs 
of complex numbers, U : H — » ffikeBzC^ . A ket localized 
at site m transforms as 



U\ma) = — ® keB z e lkl " m £ a , 

Z7T 



(Bl) 



where r m = miAj + 7712A2, and £ Q is the column vector 
of K elements with the a th entry equal to 1, and the rest 
0. 

The Fourier transform of any operator A is 



UAU* = ® keBZ A(k). 



(B2) 



The momentum space derivatives of an operator trans- 
form as 

(® keBZ d kj A{k))U=-i[r :i ,A], j = x,y, (B3) 

The following identity of traces has to hold for a clean 
infinite system 



1 



/ Tr{A(k)}rf 2 k = lim ]tTt a {A}, (B4) 
J A-yoo J\ 



(2n)- 

where Tr is the trace over orbitals, whereas Tr^ is the 
trace over orbitals and over the sites included in a patch 
of the system of area A. 

Spatial disorder configurations 5 are distributed ac- 
cording to a probability measure dfi(6), which by as- 
sumption obeys the following properties. The probabil- 
ity measure for disorder configurations is invariant under 
spatial translations (homogeneity); any subset of the dis- 
order configuration space invariant under translations is 
a set of measure 0. These properties amount to requiring 
that the probability measure for disorder configurations 
is ergodic with respect to spatial translations. Ergodic- 
ity implies that spatial and disorder configuration aver- 
ages are interchangeable. The averaging in Eq. ( 53 ) can 



be performed over a small system with many disorder 
configurations to the same effect as on a system in the 
thermodynamic limit with a single disorder configuration 
(self-averaging property). 

For each disorder configuration 6, one defines an oper- 
ator Ag (for example, the Hamiltonian Ms). One cannot 
define k-space calculus rules because the Fourier trans- 
form is no longer defined, but one may replace calculus 
rules in the Brillouin zone by the "non-commutative" rule 



and integration becomes 

lim ]rTv A {A} = f dn(5)tr {A s }, (B6) 

where the trace at the origin is given by tro{A} = 
Tr{7To^47To} involving the projector ttq = J2 a |0,a;)(0,a| 
onto quantum states at the origin m = (0,0). 

With the help of such calculus rules, the Chern number 
in the disordered continuum becomes 



v{E)=2m J dn(5)tT {P g (E)[-i[. 



;,P s (E)},-i{y,P s (E)}}}. 

(B7) 

Ref. |101j contains a proof of the fact that this expres- 
sion for an infinite system is the analytical index of a 
Fredholm operator, and thus an integer. This integer 
can only change if the energy E crosses a region of ex- 
tended states in the spectrum, i.e. something like a bulk 
band. For our numerical calculations, used the develop- 
ments in Ref. [102J, which gives a fast converging formula 
for the Chern number of Eq. (B7). The Brillouin torus 



is discretized into N x x N y patches of momentum space 
area A x A y , and integration is approximated by Riemann 
summation. 

Derivatives in k-space are approximated by a finite- 
differences formula 

Q 

dk,Pk -> S ki Pk n = ^ Cj[i\ n+J -A s - Pk a -jAi\, i = x,y 
j=i 

(B8) 

where the coefficients c m are chosen to insure exponential 
convergence in the limit of large N x>y , and Q is an ad- 
justable number of order N x<y . Then the Chern number 
becomes 

k„ 

= -27riTi{® kneBZ (Pk n [5 kl P kn ,S k2 P kn ])}.(B9) 

We can now replace the derivatives in momen- 
tum space by their real-space counterparts through 
U(e- lr -*P(E)e ir -*)W = ® k eBZ-Pk+A(£). We arrive at 
the final formula for the Chern number 



v s [E) = 
2ni 



N X Ny 



(ma\Ps{E) H Pi (E)\ ,-i[y, P s (E)\ } \ma) , 

m.(\ 

(BIO) 



where the sum is now manifestly over a finite system. 
The bracket used above is derived from the original com- 
mutator [x,Pg(E)] 

[x,P s (E)\ = 
Q 

iJ2 Cj (e- ijx - A *P 5 {E)e ijx - A °> - e ijx - A °>P s (E)e- ijx - A °>), 



3=1 



d ki A s 



-i\rj,A s ], V<5 



(B5) 



(Bll) 
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and an analogous equation holds for y. The result needs 
to be averaged over disorder configurations, 

v{E)= [dn(6)v s (E). (B12) 



Appendix C: Interacting fermions on the Kagome 
lattice 



Then at mean-field level the nearest-neighbor interaction 
amounts to shifting the bands by an amount constant in 
k-space. For strong enough interactions V, the bands 
will close gaps, but as long as V < \t\, which is the size 
of the band gap at <f> = ? , the picture cannot be changed 
qualitatively and the topological phase is robust to this 
interaction effect. The cancelation of any k-dependent 
shifts is a feature of the C3 symmetry of the lattice. 



In this appendix we analyse a Hamiltonian describing 
a fermionic nearest-neighbor interaction on the Kagome 



lattice: 



(CI) 



The Fourier transformed interaction Hamiltonian reads 



V 

N 

cos 

COS 



Ink 



' 2 
Ai 



a Ll-k a Al a C,n+k a C,n 



a C,l-k a C,iaB,n+k a B,n 



, (C2) 



where N is the number of unit cells on the lattice. We 
begin with the model containing a flat middle band at 
4> = |. Mean-field theory is performed easily if wc 

switch to the band basis (<4 k , a^, a^), where the op- 
erators create particles of given momentum in each of 
the three bands, "lower" , "middle" , and "upper" . In the 
band basis, the tight binding Hamiltonian of Eq. ([3| 
describing free particles is diagonal, 



diag(Ei(k),E m (k),E u {k)). 



(C3) 



The columns of T are the Bloch eigenvectors of Jt^. Then 
the two bases are related by the following unitary trans- 
formation, 



where 




/ T Al (k) T Am {k) T Au (k) 
T k = T Bl (k) T Bm (k) T Bu (k) 
\T C i(k) T Cm (k) T Cu (k) 



(C4) 



At 1/3 filling, or when the lower band is completely oc- 
cupied, the resulting mean-field Hamiltonian amounts to 
simply 



jfMF _ a] k aik + const, 



(C5) 



keBZ 



while at 2/3 filling, when both the lowest and the flat 
band are filled, the mean-field Hamiltonian takes the 
form 



keBZ 



(C6) 



Appendix D: Bogoliubov transformation 

In this appendix we present a way to diagonalizc the 
Hamiltonian of Eq. ( 60 1 . The Hamiltonian matrix can 



be written in the following basis 

* k = (bAk&Lk&Bk&Lk&Ck^ 



(Dl) 



with the 6 Qk operators as defined on each sublattice a £ 
{A,B,C} in Eq. ((59 1. 



= £ * k ^k$k, 



keBZ 

where now the Hamiltonian may be expressed as 

/ h AA {k) h AB (k) h AC (k) 
J#tk = h BA (k) h BB (k) h BC (k) 
\hcA(k) h CB (k) h cc {k) 

whose elements are 2x2 matrices defined as 





(D2) 



(D3) 



h-aa (k) 



h a7 tp(k) 



Un 
2 

U rip ^yaa _j_ Unp 



Pa 



-k 









a/3 



(D4) 



where is the free particle Hamiltonian in Eq. ([3]). In 
the weak coupling limit U <C \t\, the Bogoliubov transfor- 
mation can be performed as follows. The b operators can 
be transformed to quasiparticle operators b via a canon- 
ical transformation 



k^k, 



(D5) 



with <f> defined as in Eq. ( Dl ) with b replacing b where the 



condition for _B k to preserve the bosonic commutation re- 
lations between the b operators, [6 Q k, bg k ,] = <5a/3<5kk' and 
[6„k, fyak'] = is the following pseudo-unitarity condition 

B^Bl = S, where £ = id 3X 3 ® 0-3, (D6) 

where id3 X 3 is the identity acting on sublattice space and 
CT3 is the third Pauli matrix. _B k must diagonalize the 
Bogoliubov Hamiltonian of Eq. (D3); let 

(D7) 



Blm k B k 



jr k . 



From the condition in Eq. (|D6|) one can reexpress this as 

(D8) 



ES k EE^ k S k 



'k^^^ik^k — -"k ■^•'Hk&k — ^^k 

whence the matrices £J^ and are similar and it 

suffices to diagonalize the matrix Y,Jf? t to determine the 
spectrum of Bogoliubov quasiparticles. 
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Appendix E: Effective photon Hamiltonian 



state 



In this appendix we derive an effective low-energy 
Hamiltonian J^ e ^ for the photons in the Jaynes- 
Cummings lattice model deep in the superfiuid (delocal- 
ized) phase. 

The low energy state is the condensate state of n pho- 
tons and all of the 3JV two-level systems are in their 
ground state. 

\G n ) = \g n ) ® (g) U) m (El) 

rn 

where the operator a liL creates a photon in the lowest 

1 

n(Co)"io> (E2) 



,t 

band at wavevector k and 
1.9 



and |0) is the photon vacuum. 

Consider the regime of weak coupling between photons 
and two-level systems described in the main text 



g < \t\ < hio. 



(E3) 



Under realistic experimental conditions, hw rj e, which 
causes the two-level systems to be ordered in the ground 
state, and motivates our approach and to determine an 
effective photon-photon interaction starting with the ket 
in Eq. (El). Defining the projector on the ground state 



and the projector onto all excited states of the system as 



Q=\G„)(G n \ and If = id - Q. 



(E4) 



These operators project into the subspace with n polari- 



aLa 



a m a m 15 the polariton number. 



tons, where 

The indices m here represent the lattice coordinates. 
Then 



J%Q = nEi(k = 0)Q + 

^2ga m \g,i)(g n \ ® | t)mU 

m 

QJ% = nEi(k = 0)Q + 

^g\gn){g n W m ® | |) m (t 



4-)n(-i- |n j 



-l-)n(4- In- 



(E5) 



Above, Ei (k) is the dispersion of the lowest band. Simi- 
larly we will denote by E m ^ u the dispersion relations for 
the middle and upper bands, respectively, of the original 
tight-binding Hamiltonian. The following operators are 
now completely specified using the equations above 



QJ%Q = nEi(k = 0)Q 

KJ%Q = J%Q - QJftQ 
QJ%tK — QJ%t — QJ%IQ 



(E6) 



Finally the projections onto excited states can be ex- 
pressed simply in terms of the projector onto the ground 



■fflp = X! ^m.naJnGn ~ JlJSj (k = 0)Q 
m,n 

III 

^9 = 9 CT m a m + CT~ 
m 

-Q,M - Jf t Q + 2QJf t Q. (E7) 
The effective low energy Hamiltonian is written as 

1 



Jf e ff = Qj%Q + QJ%K 



E-KMiK 



KJ? t Q. (E8) 



A perturbation series in the small parameter g can be 
obtained by expanding the denominator: 



1 



1 



E-KJ^ t K E - (3V p + ,%? e + Jt° t 



e-(jf p + jr e ) id 



id 



E-(J? P + J%) 



E - {,J4° p + ,J4° e ) 
The main step is to invert 



0(g 3 ) 



where 



E-(jf p + JV e ) 



= X! et7 m°r 



(E9) 



(E10) 



(Ell) 



There is a basis in which the operator J^ v + J4? e is diag- 
onal: 

a| k a ;k _E ; (k) + a f mk a mk B m (k) + a| ik a uk £ , „(k) 



-n£ i (k = 0)Q + £ ^ ( T> n 



(E12) 



where we have introduced the operators creating a parti- 



cle in the lower, middle and upper band cLS tt^j^j d Li, and 
a «k' respectively. These band basis operators are related 
to the original operators a^ k etc. by a linear transforma- 
tion which is written explicitly in Appendix [C] The cor- 
responding band dispersion relations are E\ TO)U (k), and 
this form of the tight-binding Hamiltonian is simply the 
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diagonal form of Eq. ([3| . We have now obtained a sum 
of decoupled Hamiltonians; the tight-binding part is al- 
ready diagonal in the band basis introduced here; the 
two-level system piece is already diagonal in coordinate 
space. 

Let us compute the effective Hamiltonian up to lowest 
order contribution in g, from Eq. (E9): 



E 



1 



KJ%Q. (E13) 



Only the terms containing azk=o contribute due to the 
projectors Q at the two ends of the second term. Then we 
only need to take the following terms for the evaluation 
ofEq. pl3|: 



3V f 



- a\ ai Ei(0) - e^a^ 



(E14) 



We have cast the denominator into a form which is di- 
agonal with respect to both photon quantum numbers 
and two-level system quantum numbers. The | |)(t | of 
Q,3rf?K and the | t) (I I °f KJVQ acting on the diagonal 
matrix in between them will yield the projector onto the 
down state | I) (I \, i.e. the second term in Eq. (E13) 
will be a c- number times Q. More explicitly, 



QjfK 



1 



E - a\ Q a l0 Ei(G) + nE^Q - e£, 



Om<7 m 



KJ^Q 



\9n){9n 



9 2 E 



m flm a m 



E - 4,0,0^(0) 

,2 



\9n)(9n\®6§\l)m(l 



ng- 



E-{n- 1)^(0) 
Then to lowest order g 2 

M' eS =Q (ei(0) 



Q 



9 2 a\ m 



Hamiltonian: 



3#> e ff= [Ei(0) + 



9_ 

Ei(0) 



(E17) 



which shows that at lowest order the interaction shifts 

2 

band minimum by an amount El ^_ e ■ 

Higher order contributions yield effective interactions 
for the photons. The next nonvanishing contribution is 
in g 4 , and comes from the term 



Jf eff D QJ%K 



1 j // 

E-(je p + je e )J ' 9 



i 



E-<& + jQj je > K * lQ (E18) 
After a similar calculation to the one above, we obtain 
the following form for this term, 



D 



9_ 

(2E l (0) 



E 



t t 



(E19) 



In the main text, we have denoted the ground state en- 
ergy of the single particle tight-binding model as Ei (0) = 
E . Peeling off the projector onto the ground state, Q, 
we may rewrite the effective Hamiltonian Jt?^ as 



E a 



E 



ho 1 t 



kg" 



(2E - c )3 



al a w al a l0 , (E20) 



E- (n- 1)^,(0) 



where /a and f± are dimensionless quantities defined as 
-Q. (E15) f 2 = (\T Al0 \ 2 + |T s;o | 2 + |T c;o | 2 ) and / 4 = (|T A/0 | 4 + 
|Ts;o| 4 + |7c7o| 4 ), where the T's are the elements of 
the unitary transformation from sublattice basis to band 
basis, introduced in Appendix [C] a Q k = T a ikaik + 



E-(n-l)E t (0) 



(E16) 

To retrieve the effective Hamiltonian for the theory, we 
need to take the limit E — > nEi(k — 0), i.e. this is 
only a valid expression close to the state ground state 
onto which Q projects. Taking the limit above gives us 
an effective shift of the k = energy. Peeling off the 
projector Q, we obtain the following effective low-energy 



TamkaLk + T aukal lk for each a = A, B, C. 

Under realistic experimental conditions, where both Eq 
and e are in the microwave range, one can obtain a repul- 
sive interaction with g 4 /(2Eo — e) 3 > 0. Since we have 
only considered the ground state in which all photons are 
at wavevector k = 0, this approach does not provide pre- 
cise information about the range of the interaction, but 
it is telling that a Hubbard-type repulsive photon-photon 
interaction can emerge. 
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